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ABSTRACT 

This report documents the procedures and results of a study of a conceptual system for measuring 
the debris environment on the Space Station. The study was conducted in two phases: the first 
consisted of experiments aimed at evaluating location of impact through panel response data col- 
lected from acoustic emission sensors; the second analyzed the available statistical description of 
the environment to determine the probability of the measurement system producing useful data, and 
analyzed the results of the previous tests to evaluate the accuracy of location and the feasibility of 
extracting impactor characteristics from the panel response. The conclusions were that for one 
panel the system would not be exposed to any event, but that the entire Logistics Module would 
provide a modest amount of data. The use of sensors with higher sensitivity other than those used 
in the tests could be advantageous. The impact location could be found with sufficient accuracy 
from panel response data. The waveforms of the response were shown to contain information on 
the impact characteristics, but the data set did not span a sufficient range of the variables necessary 
to evaluate the feasibility of extracting the information. 
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1.0 INTRODUCTION 

While developing techniques for repairing damage to the Space Station by meteoroid or de- 
bris impact, a need arose for a system to locate perforations that may be small and hidden. 

A previous NASA proposal to develop a "nervous system" for the Space Station (Ref 1), 
included the suggestion that acoustic emission sensors could be used to locate impacts. It was 
stated that acoustic emission signatures might be correlated to specific types of event, but that sen- 
sor development was inadequate to allow the extraction of impactor characteristics. 

A preliminary experiment conducted at BAC in 1986 on IR&D funding showed that acoustic 
emission sensors would provide a practicable means of collecting data from the impact response of 
a panel, and a simple algorithm was able to determine the impact location with satisfactory pre- 
cision. This led to the present effort aimed at evaluating the potential of an impact detection and 
characterization system added on to an existing NASA/MSFC contract in two phases. 

The purpose of the first phase was to evaluate the limitations of the location method from the 
point of view of designing a system with a reasonable number of sensors. The effort consisted 
mainly of conducting experiments. 

The purpose of the second phase was to develop a conceptual design of an instrumentation 
system for deployment on the initial Space Station Modules for measuring, mapping, and charac- 
terizing the meteoroid/debris space environment. The data from this instrumentation would be used 
for the design of the "evolutionary" Space Station module meteoroid/debris protection system. 

This report summarizes work accomplished under the NASA Contract NAS8-36426, 
"Integrated Wall Design and Penetration Damage Control", during the periods 1 January through 
May 31 and 1 August through 16 December 1987. This program was sponsored by NASA 
Marshall Space Flight Center with Mr. Sherman Avans as the Contracting Officer's representative. 
Performance of this contract was under the direction of the Engineering Technology Organization 
of the Boeing Aerospace Company with Mr. Paul Stem as Program Manager and Mr. Alex 
Coronado as Technical Leader. 

A report (Ref 2) released at the completion of the first phase is superseded by this one. 
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2.0 THE MEASUREMENT SYSTEM 

A concept for a measurement system was drawn. Data from the existing data-base on space 
debris distributions were used to identify the probability of the concept acquiring meaningful im- 
pact data. Data from impact experiments were then used to establish the feasibility of the system. 

Two aspects of the proposed measurement system were evaluated: the sensitivity of the sen- 
sors in relation to the probabilities of experiencing measurable impacts and the ability to extract the 
needed data from sensor measurements. For these evaluations, the existing data base on space 
debris distributions were reviewed, and impact experiments were performed on panels, the mo- 
tions being measured with several sensor types. 

2.1 The Measurement System Concept 

The concept for a space experiment uses an instrumented panel as a particle impact detector. 
Impacts of space debris on a panel will produce a vibrational response which can be monitored by 
sensitive yet simple instruments. Since the motions propagate as waves, relatively few sensors are 
needed to monitor a large area, typically one per square foot or more. Thin film sensors, cheap and 
easy to use and repair, can be placed on the inside of the shield plates, with electrical connections 
to a simple data collection system. The data extraction methods discussed below derive from 
waveform analysis and are not dependent on amplitude calibration, even if this changes with 
temperature. 

By using existing panels such as the meteoroid/debris shields, the experiment is cost effective 
and not intrusive. The panels can be recovered eventually for analysis of any impact craters. The 
shield panels on the Logistics Module of the Space Station provide a useful opportunity for this 
experiment: they will be shielded from meteoroids from above by the rest of the structure, the 
surface is cylindrical in the plane of orbit and thus provides a directionally sensitive experiment, 
and the panels will be recoverable every 6 months when the module is retrieved. They are ap- 
proximately 55 in. (1.4m) square with an area of about 2 sq.m, each, as sketched in Figure 1. 

The data collection system is described by the electronic schematic given in Figure 2 which 
shows the functional blocks needed to provide the basic data functions of the system. These begin 
with a collection of sensors, typically 10 to a panel on each of up to 10 panels. For each sensor 
there is a microchip to provide digitization and a data buffer for about 200 points. Data flows 
through the buffer into a trigger circuit which selects the neighbors nearest to a given sensor which 
triggers at a preset signal level. This feature will only be necessary if required by limitations on 
computation or transmission systems. Data from these neighbors is then drawn from the buffers 
and recorded for further processing in-situ or for transmission. 
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Figure 2. Schematic of Data Flow 


4 








D 180-30708-1 


The digitized data buffers will store approximately 0.5 ms of data at a sampling period of 5 
microsecs. This is processed in real time to trigger recording at the first occurrence of a preset sig- 
nal-to-noise ratio from any of the sensors. Signals from the immediate neighbors of the trigger 
sensor will be recorded for a period of 1 ms. 

These signals will be used to locate the impact and to evaluate characteristics of the impact, eg 
its amplitude in time and spectral domains, its rise-time, and its phase spectrum. Eventual recovery 
of the panels will allow examination of the impact craters for further analysis of the impactor. 

2.2 Probability Evaluations 

The probability of impacts which could be monitored by the panel response system was esti- 
mated. To estimate the limits of detection, the system was assumed to respond to the transverse 
motions of elastic bending. These are governed by a dual wave equation whose solutions are 
determined by imparted momentum. For the small particles discussed below (with diameter-to- 
plate thickness ratio less than 0.1) the imparted momentum would be determined by the projectile 
momentum. The majority of the tests described in this report were in the elastic regime. The 
smallest impacts in the experiments used BB pellets. These were .172 (4.4 mm) dia. steel balls 
(density = 7.8 g/cc) having a mass of 342 milligrams. Velocities were around 280 ft/sec. (85 m/s), 
thus momentum was about 85 g.m/s. A set of tests was made using .156 in (4.0 mm) dia. nylon 
spheres (density = 1.1 g/cc) having a mass of 37 milligrams, with velocities about 540 f/sec (165 
m/s), and momentum of 6 g.m/s. The signals from the AE transducers were in most cases quite 
high so that it is plausible that a detection limit using these transducers would be a momentum of 
about 0.1 to 1 g.m/s. 

The prescribed Space Debris environment for the Space Station as given by D.J.Kessler 
(Ref. 3) and reproduced in Figure 3a, shows that particles up to 0.03 mm can be expected at 400 
Km and up to 0.06 mm at 500 Km altitudes. (Space Debris data presented by F. Horz (Ref 4) re- 
produced in Figure 3(b) shows that particles up to 0.3 mm dia can be expected.) The most proba- 
ble velocity is 12 to 14 Km/sec as shown in Figure 4. The momentum of an aluminum particle of 
this size and velocity is about 0.55 g.m/s (or 550 g.m/s for the Horz data). 

The expected number of impacts on the sensor during a typical mission within a range of 
momenta was determined with the following ground rules: (1) The sensor will be located on the 
Logistics Module which can be represented as a right circular cylinder whose axis lies along the 
local vertical; (2) The module is not shielded from the debris threat by any other component of the 
space station; and (3) The environment of JSC 20001 (Ref. 3) was used with a factor of 2 increase 
in flux. The increased flux level was proposed for this study by D .J. Kessler, the author of 
Reference 3, to account for the anticipated increase in the debris environment after the mid 1990’s. 
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Figure 4. Probability Distributions for Debris Velocity 
(According to Kessler (JSC 2000 1)) 
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The derivation of the governing equations is given in Appendix B. The two equations for 
the probability of momentum m or greater and the expected number of impacts with momentum m 
or greater were solved by the numerical approximations given in Figures 5 and 6 respectively. 
These equations were programmed using existing subroutines from the BUMPER analysis where 
possible. The momentum distribution is a function of the minimum diameter chosen, therefore, the 
minimum diameter was included as a user input It was also desirable to include a multiplier on the 
flux equation as a user input. The rest of the inputs relate to the sensor geometry and exposure 
time as well as the integration steps. 


G(m) = P(MZm) = j 


NPHI NT 
/*! 4*1 


m 


CV(9;)[COS(0; - <t>, )] + ^ 
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[cos(0, - ^)]V(0,.)A0A4. 


where __ 

H = 1 -H 

— _ F DEB (D ) 

" F DEB (Do) 

NPHI - Number of ()) steps 

NT = Number of 0 steps 

Figure 5. Probability of Occurrence of Momentum m or Greater 
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= Exposure time of sensor 


Figure 6. Expected Number of Impacts of Momentum m or Greater 
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The analysis was performed using the default inputs for both the normal and absolute mo- 
mentum calculations. The results are shown in Figures 7 to 9. The results indicate that for a typi- 
cal 6 month mission a sensor in the most advantageous location will encounter fewer than 1 particle 
with sufficient momentum to register. Estimates for the expected number of impacts on the entire 
Logistics Module, as a function of particle diameter are given in Appendix B. The results show 
that for a momentum measurement system with a threshold of about 0.1 g.m/s, the probability of 
an event in the six-month period of each experiment is about 0.3. 



Figure 7. Momentum Distribution, M <m 
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Figure 8. Momentum Distribution, M>m 



Figure 9. Expected Number of Impacts, M> m 
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3.0 IMPACT TESTS 

Tests were made on a 1.22 m x 1.83 m x 4.76 mm (4 ft x 6 ft x 3/16 in) thick aluminum 
panel. The impactor was either a copper clad steel or a nylon BB from a common BB gun. A 
sample of the integrally-stiffened machined Space Station waffle panel shown in Figure 10 was 
used for a brief series of tests including 4 penetration tests under high velocity impact from a pow- 
der gun. 

The first few tests were made using two wide-band acoustic emission sensors with a two- 
channel digitizer. A multi-channel device was acquired later to enable the recording of signals from 
several transducers on each test. 

Broad band acoustics emission transducers were used for most of the measurements, but thin 
film polymer sensors and a low-frequency piezoelectric sensor were also used. 

A large low-frequency (100 KHz) ultrasonic transducer was placed under the impact point in 
several tests to provide an impact time fiducial. 

3.1 Descriptions Of Tests 

The sequence of tests that were conducted to gather relevant data on sensor locations and 
gains, and impact locations and velocities is given in Tables 1-4. The test configurations of the 
later tests are shown in Figure 11a through f. The purpose of each test and a discussion of the test 
progression follows. The results are presented and discussed below. Reproductions of all recorded 
waveforms are given in Appendix A. 

The first three tests of 01/21 were made to examine changes in the waveforms between two 
sensors spaced apart in a line. They were placed at 102 mm (4 in) and at 254 mm or 381 mm (10 
in or 15 in) from the impact. Waveforms for the three tests are shown in Appendix A as Figure A- 
1. The several waveforms at 102 mm (4 in) and those at 254 mm (10 in) are quite similar to each 
other indicating repeatability. The 254 mm (10 in) waveforms exhibit the expected dispersion 
(spreading of the pulse to lower frequency at longer time) but their amplitudes were not the same. 

Repeatability between tests was of concern, since distance effects could only be evaluated 
from data on repeated tests, because only two recording channels were available 

The next series of six tests was made on 02/10 to provide data at several distances, with im- 
proved experimental techniques. The two transducers for each test were bonded with Loctite and 
the launcher was mounted at a fixed distance from the panel. The data for sensors at 125 mm, 178 
mm, and 254 mm (5 in, 7 in, and 10 in)(figure A-2) showed expected changes with distance. The 
data from the near-in reference sensor at 76 mm (3 in) (figure A-3) showed considerable variation 
in signal strengths and waveforms among the tests. 
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Complete Waffle Panel 



Closeup of shot (05/05) 

Figure 10. Waffle Panel After Hypervelocity Impact 
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Table 1. Preliminary Tests 


Note: Digitization Interval: 0.5 Microsec 


Date 

Data Id 

Dispersion ' 

Data 

01/21/87 

IMP1.DAT 

01/21/87 

IMP2.DAT 

01/21/87 

IMP3.DAT 

02/10/87 

IMPACT.DAT 

02/10/87 

IMPACT.DAT 

02/10/87 

IMPACT.DAT 

02/10/87 

IMPACT.DAT 

02/10/87 

IMPACT.DAT 

02/10/87 

IMPACT.DAT 

Impact Velocity Data 

03/05/87 

0306.AE 

03/06/87 

0306.AE 

03/06/87 

0306.AE 


Evaluate Transducer #4 
03/20/87 0320.AE 
03/20/87 0320.AE 


a - Gain = 2.0 


Sensor Distance 
From Impact (ins 


#1 #2 #3 #4 



Table 2. Multi-Channel Fixed Distance Data 


Note: Digitization Interval: 1.0 Microsec 


Date 

Data Id 

Sensor Distance 
From Impact (ins) 

Measured Impact 
Velocity 
(ft/sec) 

#1 

#2 

#3 

#4 

#5 

03/25/87 

0325_1.AE 

12 

12 

12 

12 


270.0 

03/25/87 

0325_2.AE 

12 

12 

12 

12 


268.7 

04/01/87 

0401.AE 

12 

12 

12 

12 

- 

265.5 

04/01/87 

0402.AE 

12 

12 

12 

12 

- 

261.6 







0325_1.AE 

1.56 

1.56 

1.56 





0325_2.AE 

1.56 

3.13 


3.13 

1.56 



0401. AE 

1.56 

1.56 

3.13 


- 



0402.AE 

3.13 

3.13 


3.13 

- 
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Table 3(a) Location Tests 
Notes: Digitization interval: 1.0 microsec 

Locations of sensors and impacts shown in Figs 11(a) to (f). 


Date 

Data Id 

Measured Impact 
Velocity 
(ft/sec) 

Impactor 

Flat Panel 

04/06/87 

0406_1.AE 

260.3 

Steel ball 

04/06/87 

0406.2.AE 

260.3 

Steel ball 

04/06/87 

0406_3.AE 

270.2 

Steel ball 

04/06/87 

0407- 1.AE 

268.7 

Steel ball 

04/14/87 

0414_1.AE 

520.5 

Nylon sphere 

04/14/87 

0414.2.AE 

540.8 

Nylon sphere 

04/14/87 

0414.3.AE 

520.2 

Nylon sphere 

04/14/87 

0414_4.AE 

526.0 

Nylon sphere 

Waffle Panel 
Low velocity 

04/20/87 1 0420_1.AE 1 

536.0 

Nylon sphere 

04/20/87 

0420_2.AE 

539.7 

Nylon sphere 

04/20/87 

0420_3.AE 

537.8 

Nylon sphere 

04/20/87 

0420_4.AE 

538.9 

Nylon sphere 

04/20/87 

0420.5.AE 

555.1 

Nylon sphere 

04/20/87 

0420.6.AE 

520.4 

Nylon sphere 

04/23/87 

0423_1.AE 

262.1 

Steel ball 

04/23/87 

0423_2.AE 

263.4 

Steel ball 

04/23/87 

0423_3.AE 

256.1 

Steel ball 

04/23/87 

0423_4.AE 

267.9 

Steel ball 

04/23/87 

0423.5.AE 

276.9 

Steel ball 

04/23/87 

0423_6.AE 

258.4 

Steel ball 

High velocity 

05/04/87 | 0504_1.AE 1 

3.10 km/sec 

1/4" dia x 1/4" long Al, no MLI 

05/04/87 

0504_2.AE 

3.22 km/sec 

1/4" dia x 1/4" long Al, with MLI 

05/05/87 

0505_1.AE 

1.27 km/sec 

1/4" dia x 1/4" long Al, no MLI 

05/06/87 

0506_1.AE 

3.30 km/sec 

1/8" dia x 1/8" long Al, no MLI 


Note: MLI=Multi-layer insulation 
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Table 3(b). Sensor Gains and Types for Location Tests 


Test Number 

Sensor Gain (x 1 E04) 

1 

2 

3 

4 

5 

6 

7 


3.13 

3.13 

3.13 

3.13 





1.56 

7.81 

3.13 

1.56 




0406-3 

3.13 

7.81 

3.13 

3.13 




0407-1 

3.13 

7.81 

1.56 

3.13 




0414-1 

0.78 

0.78 

1.56 

0.78 




0414-2 

0.78 

0.78 

1.56 

0.78 




0414-3 

0.78 

0.78 

1.56 

0.78 




0414-4 

0.78 

0.78 

1.56 

0.78 




0420-1 

1.56 

1.56 

0.016(K) 

1.56 

0.078(K) 



0420-2 

1.56 

1.56 

0.160(K) 

1.56 

0.078(K) 



0420-3 

1.56 

0.78 

0.16-(K) 

3.12 

0.078(K) 



0420-4 

0.16 

3.12 

0.03 1(K) 

0.31 

0.016(K) 



0420-5 

0.78 

1.56 

0.07 8 (K) 

0.78 

0.03 1(K) 



0420-6 

0.78 

0.78 

0.016(K) 

0.78 

0.078(K) 



0423-1 

7.81 

3.13 

3.13 

7.81 

0.16(K) 

3.13(S) 

0.078(S) 

0423-2 

7.81 

7.81 

1.56 

7.81 

0.1 6(K) 

3.13(S) 

0.078(S) 

0423-3 

7.81 

7.81 

1.56 

7.81 

0.1 6(K) 

3.13(S) 

0.078(S) 

0423-4 

7.81 

7.81 

1.56 

7.81 

0.1 6(K) 

3.13(S) 

0.160(S) 

0423-5 

0.078 

7.81 

0.78 

0.78 

0.1 6(K) 

1.56(S) 

0.078(S) 

0423-6 

0.078 

3.13 

0.78 

3.13 

0.1 6(K) 

3.13(S) 

0.078(S) 

0504-1 

7.81 

7.81 

7.81 

7.81 

1.56(K) 

7.81 

7.81(S) 

0504-2 

15.63 

7.81 

7.81 

7.81 

1.56(K) 

15.63 

7.81(S) 

0505-1 

15.63 

7.81 

7.81 ! 

7.81 

0.78(K) 

15.63 

0.78(S) 

0506-1 

7.81 

7.81 

7.81 

3.13 

0.78(K) 

7.81 

3.13(S) 


NOTES: 

(K) designates “Kynar” thin foil sensor 

(S) designates “Sondicator” low frequency piezo sensor 

All other data are for acoustic emission (AE) sensors 
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Table 4. Impact and Transducer Locations and Edge Echo Paths 


Direct Min Echo | First Pulse 

Path (cm) I Dim 


Pathfmm 




164.92 

586.80 

722.62 

453.34 


68.07 

95.80 

119.05 

79.57 


194 

271 

297 

240 

161.2 

5.40 


680.12 

76.79 

630.72 

927.35 


90.91 

60.62 

121.88 

139.48 


567.30 

211.62 

422.02 

678.85 


97.19 

81.46 

101.45 

114.79 


263 

195 

229 

283 

152.6 

5.23 


535.49 

93.09 

677.73 

860.83 


76.91 

55.81 

117.06 

128.45 


282 

194 

305 

334 

179 

5.43 


747.41 

537.54 

112.54 

539.12 


126.17 

114.70 

68.32 

86.48 


Extrap event time 161.8 mus 


#414-4 


747.41 

537.54 

112.54 

539.12 


747.41 

537.54 

112.54 

539.12 


125.92 

114.74 

68.39 

86.11 


126.17 

114.70 

68.32 

86.48 



E 

298 

N 

265 

S 

177 

S 

259 


155.5 


5.12 

N 

414 

N 

272 

S 

183 

S 

270 


125.1 


3.03 
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(b) Tests of 04/20 and 04/23 




Figure 1 1. Schematic of Test Configurations 
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Figure 1 1. Schematic of Test Configurations (Continued) 
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Figure 1 1. Schematic of Test Configurations (Concluded) 
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The third set of three tests on 03/05 and 03/06 were instrumented for measuring the impactor 
velocity to determine whether this could be a source of variation from shot to shot. It was found to 
be quite repeatable as shown in Tables 1, 2, and 3. The waveforms at 125 mm, 178 mm, and 254 
mm (5 in, 7 in, and 10 in) (figure A-4) were similar to those of the previous tests, while those at 
76 mm (3 in) appeared to be somewhat similar. 

Three tests of 03/20 were intended to evaluate the transducer used at 254 mm (10 in) because 
doubt arose as to whether its response was inverted. The onset of the high frequency perturbation 
in relation to the low frequency appeared to be inconsistent with that of the other tests. This was 
subsequently attributed to the arrival of an echo from the edge of the plate at about the same time 
for all transducers. Two transducers were placed side-by-side at about 305 mm (12 in.). The re- 
sults (figure A-5) showed only moderate differences between them, but both sensors exhibited 
waveforms that were less characteristic of dispersive bending than in the previous tests. It is 
thought that the mass of the two together inhibited bending. It could not be clearly established 
whether there was a sign reversal. 

A preliminary test was made on 03/25 using a multi-channel digitizer to compare the signals 
from four transducers at one distance. Also, a projectile timing device consisting of two screens 
of fine wire was placed between the gun barrel and the panel to determine whether the gun perfor- 
mance was repeatable. The transducers were placed next to each other on a line perpendicular to a 
radius from the impact. The signals (figure A-6) were quite similar to each other in the first 100 
microseconds, or so, but became dissimilar after that. These waveforms exhibited little bending 
dispersion, again probably due to the mass of four transducers close together. 

Several tests were made on 3/25, 4/1, 4/2, 4/6, and 4/7 using a configuration with four 
transducers placed at 45 deg. intervals around an arc of 305 mm (12 in) radius centered on the 
impact point. An ultrasonic transducer was used at the impact point to provide a trigger in the first 
two tests (though its signal was not recorded) but it was not installed for the other tests, so 
triggering was taken from the other transducers. The data are shown in figures A-7 to A- 13. The 
first two of these showed that transducer #4 was still suspect, as its signals were low by a factor of 
10. Indeed it was found to contain a faulty connector. Problems were encountered with the 
projectile timing device and were resolved during this sequence of tests. The transducer closest to 
the impact (no. 2) exhibited unexplained high spikes at an early time in the last three tests. 

Four tests were then conducted on 04/14 to obtain data similar to those of the original IR&D 
effort, but now using multi-channel recording of four transducers arranged in a 305 mm (12 in) 
square, so that all data were simultaneous. The projectile was a 4.0 mm (5/32 in) dia. nylon 
sphere. The data given in figures A- 14 and A- 15 show lower signals than for the previous tests 
with the copper-clad steel BB, but apparently with a higher content of high frequency. 
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Six tests were made on 4/20. Additional transducers were used to investigate alternatives to 
the acoustic emission sensors. These included a low frequency bending sensor ("Sondicator") 
based on a 25 KHz crystal, and a thin piezoelectric film ("Kynar"). Data are shown in figures A- 
16 to A-21. 

Tests on the waffle panel were begun with six low velocity (BB) shots on 4/23. The sen- 
sors were arrayed within the waffle squares to provide transmission directly within a panel as well 
as across the stiffener ribs. Data are given in figures A- 22 to A-21. The signals varied widely 
among the transducers, though all showed some of the expected features. Many included high 
spikes which offset the automatic scaling, but the signal appears otherwise normal. The problem is 
thought to lie in the digitizer which uses an interpolation procedure to fill in data between sampling 
times. 

High velocity tests on the waffle panel were conducted on 5/4, 5/5 and 5/6 using a BAC two- 
stage powder gun with cylindrical aluminum pellets as illustrated in Figure 12. The installation of 
the shield plate and the multi-layer insulation materials is illustrated in Figure 13. These tests were 
designed to conform to predictions of penetration made in BAC IR&D studies (reference 5) as rep- 
resented by the graph of Figure 14, showing a penetration threshold for projectile diameter and 
velocity. As indicated, test conditions were selected to provide some penetrating and some non- 
penetrating shots: one with 6.4 mm (1/4 in) dia and 6.4 mm (1/4 in) long pellet at 3 mm/ 
microsecond onto a shield of 0.063 in aluminum; the same with 20 layers of insulation on the 
panel; the same at 1.5 mm/microsecond; and one with a 3.2 mm (1/8 in) dia. and 3.2 mm (1/8 in) 
long pellet at 3 mm/microsecond. 

The first test was unsuccessful because the trigger was incorrectly timed. Further, no pene- 
tration occurred when it was expected suggesting that the fault lay in the gun firing. The remaining 
four shots were satisfactory. Signals for these shots are shown in figure A-28 to A-31. 
Photographs of the impact areas are shown in Figure 10. 

Any possible effects of air-borne shock from the powder-gun were mitigated by the shield 
plate. Induced motions would have to propagate along the shield plate, down the support screws, 
and along the panel to the sensors. The propagation would attenuate and delay the motions 
considerably so that they would not interfere with the motions from penetration. 
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Figure 12. Waffle Pane, Impact Tee, Setup for Impact Detection/Location System 
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Figure 13. Installation of Shield and Insulation on Test Panel 



Figure 14. Selected Test Points Versus Pretest Predictions 
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4.0 THEORY 


4.1 Plate Response 

The dynamic response of a plate takes the form of in-plane and lateral motions. These 
generally occur in a wide range of modes, some propagating as waves with frequency-dependent 
wavespeeds and others as stationary oscillations which do not propagate. Up to moderate frequen- 
cies the in-plane response propagates by simple waves having a constant wavespeed known as the 
plate speed. This is typically slightly less than the dilatational (longitudinal) speed for the material, 
because the plane stress state caused by lateral relief of stress at the faces of the plate is less stiff 
than the plane strain state of a plane wave. The transverse response behaves like ordinary bend- 
ing at low frequency with a diffusion-like character having an apparent 1 speed that varies with the 
square-root of frequency. At high frequency it behaves like independent Rayleigh surface waves 
on the two faces of the plate. The relative magnitudes of these various motions are dependent on 
the loading and the boundary conditions. 

4.2 The General Solution 

Solutions of the differential equations for in-plane motions, u(x,z,t) and lateral motions, 
w(x,z,t) of a thin plate are found (reference 6a) in the propagating exponential forms 

u = U(z)e ik ( x - Ct ) 


and 


w = W(z)e&( x -c0 

These propagate at the wavespeed c and oscillate in time at a frequency © = kc with spatial 
wavenumber k. The thickness variations U(z) and W(z) are determined by two thickness scale 
factors 

q = k(l-c^/ci^)l/2 

and 

s = k(l-c 2 /c2^) 1 ^ 2 


1 A motion with such wave speed dependence is not a true wave because the speed of low 
frequencies approaches zero, and no wave packets of similar frequency can group together as all 
frequencies have different speed. 


FA3£ Bum WOT FKJWcD 
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with ci and C2 being the longitudinal and shear speeds of the material, respectively. 

These quantities are related by the two equations for: 

1) in-plane motion: 

tanh(sh)/tanh(qh) = 4sqk2/(k2+s2) 
and 

2) transverse motion: 

tanh(qh)/tanh(sh) = 4sqk2/(k2+s2) 

These are referred to as the dispersion relationships which determine the wavespeed c 
through q and s as functions of k and thus of 0). Graphical representations are available as illus- 
trated in figure 15 (from reference 6a). 

4.3 In-Plane Motions 

Up to moderate frequencies, the in-plane motions are governed by the common wave equa- 
tions, in which the wavespeed is the plate velocity 

Cpl = (E / p (1 - 1)2}} 1/2 = {(1 - 2t))l/2/(l - u)} cl 

where E is Young’s modulus, p is the density, x> is Poisson's ratio, and cj_ is the longitudinal 
wavespeed for plane waves in unbounded media. For '0=1/3, typical of aluminum, the factor is 
0.904. Since cl in aluminum is about 6.3 mm/microsecond, then Cpl = 5.7 mm/microsecond. 

At high frequencies where motions through the thickness become constrained by lateral iner- 
tia effects, the waves propagate slower, at the Rayleigh velocity, which for aluminum is cr = 3.1 
mm/microsecond. 

From a point source, these waves propagate cylindrically, decaying, because of increased 
area, as the logarithm of distance. Typically this reduces amplitude to negligible levels within 20 
to 100 source diameters, about 120 to 600 mm (5 to 25 in) for a 6 mm (1/4 in) impactor. 

4.4 Lateral (Bending) Motions 

A simplified analysis based on the Timoshenko Beam Theory (reference 6b) leads to the 
approximate but explicit dispersion relationship: 

k = (to/ V2cpl){ 1 + a ± [(l-a)2 + (2cpi/r co)2] 1/2 } 1/2 
where a is the ratio of Young's modulus to a factored shear modulus, E/fG, (conveniently taken to 
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Dispersion curves for the Timoshenko theory and for the 
two lowest transverse modes in a slab and in a circular rod 
(Poisson's ratio v - 0.29) 

Figure 15. Calculated Dispersion Curves for a Plate 



Figure 16. Calculated Dispersion Curves for Bending of Aluminum Plates of Two Thicknesses 
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be the square root of the ratio of longitudinal to Rayleigh wavespeeds), and r is the radius of gyra- 
tion of the plate, r = h/ Vl2 with h the thickness. This has two modes, one for the positive and one 
for the negative radical sign. At low frequencies, the positive mode gives imaginary wavenum- 
bers, so that these modes do not propagate. Above a certain cut-off frequency, however the 
wavenumbers are real and a high speed propagation develops. 

From this formula the wavespeed c = to/k can be written as 

c= V2c p i/{l+(cpi/cR)2 ± {[1 -(c p ]/cr)2]2 + 4(cpi/c b 4 } 1/2 } 1/2 

where cb is the bending speed (Eh3/12 (J.) 1 / 4 co 1 / 2 ’ an d p- pb is the plate areal mass density. 
The coefficient in this formula is equivalent to the vibrating beam stiffness coefficient (EI/p.). The 
dispersion formula thus has three parameters, Cpi, cr (or their ratio which depends only on 
Poisson’s ratio), and h. 

At low frequencies, the negative mode reduces to c = cb, which is 0 at 0 frequency, and at 
high frequencies to c = cr. The positive mode is usually of small magnitude and so is not con- 
sidered further here, though it should not be dismissed. 

A graph showing the dispersive wavespeed for 3.2 mm (1/8 in) and 4.8 mm (3/16 in) thick 
aluminum plates is given in Figure 16. 

4.5 Transient Motions 

The motions controlled by these dispersion relations can be expressed by Fourier Integrals as 
superpositions of all frequency components. Thus, the response at the impact point can be written 
as 

u(t) = J U(co) exp(i ci)t)dco 

where U(co) is the transform of the motion, a complex variable, and can be written as 
U(co) = A(co) exp[i0 o (co)] 

A is the amplitude of each frequency component, and 0 O is the phase. For an impulsive source, 
all components start as maxima at the same phase 7t/2. When these components propagate, each 
at a different speed, the signal is given by the Fourier Integral with shifted time t' = t - x/c(co): 


I 


! 
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u(x,t) = J U(co) exp[i co(t-x/c)]dco 
so that the apparent phase of the transform of u(x,t) is 
6(G)) = k /2 - x/c(to) 

Since the time of impact is not generally known, an unknown time-shift through to occurs in 
the data, resulting in a spectral change. The measured waveform then has the spectral phase 

0ni = 0 + co to = 71/2 - cox/c + CD to 

This formula shows that the phase spectrum contains the distance and time shift explicitly. 
Since the frequency dependence of the terms is known, it should be possible to use a regression 
analysis on measured spectra to obtain these. However, the phase varies as G)l/ 2 (since the speed 
c varies as co^ 2 and this leads to an infinite slope at zero frequency. Thus, the phase spectrum 
changes fast near zero. For a numerically-derived spectrum at discrete frequencies, the change 
between each point can be larger than 2tc. 

Numerical analysis can only evaluate the principal phase between say k and -it so that uncertainty 
arises at the low frequency end of the spectrum. 

4.6 Determining Location 

Measurements of motion at several places provides data on time of propagation at a known 
speed from the unknown impact point to several known places, but with an unknown time shift 
common to all sensors. There are three unknowns: the time of impact, and the two coordinates of 
the impact. Two methods are available for determining the unknowns: a minimization of the 
squared error, and the intersection of difference hyperbolae. The former is useful for many sen- 
sors, and the latter for few. 

a.) Minimized Error Method 

For a known point (x,y) the travel time relative to an arbitrary time with the unknown shift 
to is t. For a wave travelling with a speed c from the impact point (xo,yo)> these are related by 

r = {(x-x 0 ) 2 +(y-y 0 ) 2 } 1/2 = c(t+to) 
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Any error in this determination can be written as 
e = r - c(t+to) 

and a least squares definition of best fit defined as 
Ze^ = Min 

The variables in this minimization are the unknown impact position (xo,y 0 ) and the time shift 
to- The minimization leads, by differentiation, to the equations 

Xe0e/9x o ) = Xe(9e/9y 0 ) = Xe(9e/9to) = 0 

which can be written as the following implicit equations 

x 0 = {Zex/r}/{Ie/r} 

yo = {£ey/r}/{Xe/r} 

to= Ir/c-Xt 

An iteration scheme based on these equations serves to determine xq, yo. and to from sev- 
eral waveforms, 

(b) Difference Hyperbola Method 

Since the unknown impact time is common to all sensors, it can be eliminated by considering 
the differences in time of arrivals between two sensors. Consider the time of arrival at sensors i 
andj: 


Aij = ti - tj 

The locus of all points for which the difference in distances to the two sensors is the hyper- 
bola: 


ax^ + by2 + cxy + dx + ey + f= 0 
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with 


a = (xi-xj) 2 - Aij 2 
b = (yi - yj ) 2 - Aij 2 
c = (xi - xj) (yi - yj) 

d = (xi - xj) (xi 2 - xj 2 + yi 2 - yj 2 ) - Aij 2 (xi - xj) 
e * (yi - yj) (xi 2 - xj 2 + yi 2 - yj 2 ) - Aij 2 (yi - yj) 
f = (xi 2 - Xj 2 + yi 2 - yj 2 - Aij 2 ) 2 - Aij 2 (xi 2 - yi 2 ) 


Two such hyperbolae, constructed from two pairs of sensors, will define a unique intersec- 
tion point from which the distances to the four sensors satisfy the observed differences. It is thus 
a matter of solving simultaneously two such hyperbolic equations. 
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l 

j 5.0 TEST RESULTS 

i The recorded waveforms have been analyzed to evaluate four factors: 1) the accuracy of 

i defining time-of-arrival of the signals; 2) the accuracy of impact position estimates; 3) the relation- 

ship of the Fourier transform to impact characteristics; and 4) feature analysis as a means of de- 
terming impact characteristics. 

5.1 Time Of Arrival (TOA) 

The onset of motion (or time-of-arrival (TOA), or first break in the signal) is needed to esti- 
mate wave travel time for determining impact position and for isolating the various response pulses 
for Fourier analysis. It is defined as the time at which the signal exceeds a threshold. The practical 
problem is to distinguish between signal, early or late, and noise, since the first motion is small 
and larger motions represent modes at varying wavespeeds. The fastest waves are the in-plane 
: motions which decay rapidly with distance, followed by the high frequency Rayleigh waves. 

Trailing these at ever lower speeds are the bending waves, which can be the largest 

A typical waveform (from the sensor at 254 mm (10 in), file IMP2.DAT of 1/21) is shown, 
together with early details on enlarged time scales, in Figure 17. The earliest pulse train, ampli- 
i fied at the lower left of the figure, appears to have a slight inverse dispersion, with frequency de- 

I creasing along the pulse from a period of about 3.5 microsec for a half wave (i.e. at a frequency of 

0.14 MHz) to about 1.5 microsec (0.3 MHz). The large later pulse train shown at the lower right 
1 is strongly dispersive with the frequency increasing steadily by a factor of ten or so from 3 mi- 
crosec per half wave (0.16 MHZ) to about 30 microsec (16 KHZ) after about 0.3 milliseconds. 
The clipping of this waveform does not affect these features. This behavior, which represents 
slower propagation as frequency decreases is characteristic of bending motions. The increasing 
separation of the early and late pulses with distance is illustrated in Figure 18 for traces from the 
tests of 02/10. 

Details of four different waveforms all at one distance (sensors 1 to 4 from 0325_1.AE) are 
shown on an exaggerated time scale in Figure 19 displaying measurement repeatability to within a 
j few microseconds. For wavespeeds up to 6 mm/microsecond, however, this represents a propa- 
gation distance uncertainty of up to 20 mm. 

■ Waveforms from all tests for which signals were recorded from two sensors in line with the 

impact have been analyzed for arrival time using a sequence of thresholds. For each waveform, the 
baseline of the recordings was found for the first 100 points, and then the data were shifted to this 
baseline and smoothed by forming moving averages over 7 points. The time of first occurrence of 
a signal above each of five levels in a doubling sequence of thresholds were then found. 
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Figure 17. Details of a Typical Recorded Waveform 
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The increments of distance and time between pairs of waveforms for two in-line sensors are 
given in Table 4a of Sect 3 above for the six tests of 0325-1 and -2 (excluding the waveforms 
from sensor #4 which had low signals) arc plotted against threshold level in Figure 20. They show 
that for a threshold below about 5 volts the arrival time increment is between 60 and 65 mi- 
croseconds, implying a wavespeed of 4.7 to 5.1 mm/microsecond, whereas above 7 volts, the 
times are between 105 and 120, for speeds of 2.5 to 2.9 mm/microsecond. 

For an automated system and for a more precise measure, the TO A must be determined from 
some threshold criterion independent of sensor calibration or characteristics. Several such methods 
were tested, and one based on signal-to-noise ratio as described below was selected. 

An algorithm was set up to the determine the point at which the signal exceeds the noise 
floor. The noise baseline was defined by a running sequence of least-squares fits to the transducer 
readings for times ranging from the start of recording and extending, in turn, to each subsequent 
measured data point. 

Each such fit was extrapolated to the next data point at which the signal was defined as the 
magnitude of the difference between the measured value and the baseline extrapolation. The noise 
level was taken to be the standard deviation of the fit. The signal-to-noise ratio (SNR) was then 
used to determine the arrival of the signal as the point where it exceeded a threshold. 

An example of this determination of the SNR for four records in the test of 04/01 is shown 
in Figure 21. The sensors were all on an arc at the same distance from the impact. The figure 
shows the coincidence of the first arrivals at about 205 microsecs, with an SNR of about 5, and the 
main pulses at about 255 microsecs, with SNR exceeding 30. For most cases a value of 5 for the 
SNR was found to be sensitive to the first arrival and 10 for the main pulse. The method does not 
work for low-level data from large distances where the data was essentially noise, as could be ex- 
pected. 

TOA values were determined for the first and the main pulse from each sensor in each test 
from 04/06 to 04/14 as listed in Table 4. Using the known impact positions distance-time regres- 
sion fits were made for the groups of four in each test, to allow the event time of that test to be de- 
termined as the intercept time. The data for all sensors from all tests were then corrected to these 
zero-times, and plotted together in one graph as in Figure 22. The data fall closely on one line, 
whose slope (by a regression fit) is 5.26 mm/microsec for the measured speed of the first wave. 

5.2 Position Location Estimates 

Using an iterative algorithm that minimizes the squared error in distance-time estimates, the 
time-of-anival of 04/06 and 04/07 have been analyzed to find impact positions. The results are 
plotted in Figure 23. Accuracy of location is within about 75 mm (3 in). 
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Figure 20. Differential Time-of-Arrivai Data 



Figure 21. Signai-to-Noise Ratio, (Sensors 1 to 4 of 04/01) 
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Figure 23. Comparison of Calculated and Actual Impact Positions 
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5.3 Fourier Spectral Analysis 

Spectral analysis was made by two methods: the simple zero-crossings technique and the 
Fast Fourier Transform (FFT) technique. The zero crossing method provides a simple evaluation 
of the phase spectrum from which the requirements for the FFT can be deduced. 

5.3.1 Zero-Crossing Analysis 

The recording IMP2 of 1/21 was analyzed by first establishing a baseline and then finding 
times at which the signal crossed this baseline. These times were used to calculate frequencies and 
group velocities. 

Between each pair of zero-crossings, the signal pulse can be approximated by a sine-wave 
whose half period is the difference of the times, so that 


tp,i = 2(t x ,i - t x ,i-l) 


where tp,i is the period of the signal at the i-th pair 
t x ,i is the time of the i-th zero crossing. 

The frequency is then 


fi = l/2(t Xj i - t Xj i-i) 


Such a pulse is formed by the coincidence at that time and place of all waves near this fre- 
quency having the same phase. That is, the variation of phase with frequency is zero, which is 
the definition of a group whose velocity is the group velocity, eg = dco/dK where to is the angular 
frequency and k is the wave number, used in the solution form exp i(kx-o)x). The pulse is cen- 
tered on the transducer at a known distance x at the time l/2(t x ,i + t x ,i-l) so that the group veloc- 
ity is 


Cg = 2x/(t Xj i + t X) i-i) 


Thus, determination of the zero-crossing times allows derivation of the group velocity as a 
function of frequency. The wave number k is then obtained by the integral 

k = Jdco/cg 
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and the phase distribution is then 
0 = kx - cot 

The phase velocity is c = co/k. 

The phase spectrum calculated by this method for the main pulse of IMP2 is shown in Fig. 
24a. The phase velocity is shown in Fig. 24b; where it is compared to the theoretical dispersion 
curve. There is close agreement This method is not satisfactory for noisy data, however. 

The phase spectrum shows that the main pulse is a bending motion and that its frequency 
range is from about 0.01 MHz to 0.2 MHz. The spectrum of the first pulse spans from about 0.1 
to 0.5 MHz. 

5.3.2 Fast Fourier Transform (FFT) 

To calculate the FFT of a signal it is necessary to select the data sampling period dt and the 
duration of the signal tmax. These are determined by the maximum frequency and frequency step 
in the FFT by the formulae (Ref. 7) 

dt = 1/fmax » tmax = 1/df = 1/fmin 


For the frequency ranges given above, these are 

for the main pulse: dt = 5 microsec tmax = 100 microsec 

for the first pulse: dt = 2 microsec tmax = 10 microsec 

The frequency ranges, and hence the times, are bounds only, and wider ranges may be more 
practical. 

The superposition of multiple modes is compounded by the arrival of echoes off the edges 
of the test panel. These echoes have a similar magnitude spectrum, but with a time shifted phase. 
To minimize this problem, the time of arrival of the first echo (TFE) was determined, and all 
FFT’s were made with data for times shorter than TFE. 

To determine the TFE, impact-to-sensor distances on echo paths represented by oblique dis- 
tances to the nearest edge, were calculated and are given in Table 4 above. A wavespeed of 5.2 
mm/microsec was then used to determine the echo time of the first pulse. 

The data are further prepared for the FFT algorithm by smoothing, by applying a triangular- 
weighted moving average to the time series data and then detrending by a regression fit. This re- 
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moves the static and very low frequency components of the spectrum which arise from drift and 
not from the signal. 

The magnitude and cumulative phase spectra derived from FFT's of the first and main pulses 
of IMP2 are shown in Figure 25 and 26. The magnitude spectrum of the first pulse peaks at an 
angular frequency of about 0.7 rads/microsec i.e. at a frequency of about 0.1 MHz, lower than 
0.3 given by the zero crossing analysis. The main pulse peaks at about 0.1 rads/microsec, or 
0.02 MHz, in general agreement with the zero-crossing result. The phase spectra of the two are 
lower by a factor of two and four respectively than that of the zero crossing analysis, and are re- 
versed. Note that there is a difference in scales by a factor of 2*pi between frequency and angular 
frequency in Figures 24 and 26. Since the zero-crossing result was shown above to correspond to 
the theory, the cause of the discrepancy lies in either of two problems: the calculation of phase 
from the real and imaginary parts of the FFT; or the phase shift in the FFT caused by the uncer- 
tain time origin. 

The analysis of the phase spectrum uses an arc-tangent algorithm to derive the principal val- 
ue of phase (between -pi and +pi) from the calculated real and imaginary parts of the FFT. This 
principal value phase, as illustrated by curves for sensor No. 4 of Test 03/06 in Figure 27(a) is 
discontinuous, taking occasional jumps by about 2*pi. An algorithm is appended to remove such 
jumps, converting the phase into a continuous curve as in Fig. 27(b). The accuracy of the algo- 
rithm rests on a criterion for defining the jumps. The simplest and most effective sets a threshold 
on the difference in two adjacent phase values: theoretically such a difference should be less than 
2*pi, but in practice a value of about 4 is best. Nonetheless, as seen at high frequency for the 120 
microsec curve, spurious discontinuties can be introduced. For very noisy phase data, the proce- 
dure can give in erroneous results. 

The second factor (time shift) subtracts from the phase a term linear in frequency and pro- 
portional to the time shift: 


A0 = -coAt 


At an angular frequency of 1 rad/microsec, a time shift of 30 microsecs produces a phase 
change of 30 rads, which is of the order of the discrepancy. The curves of Fig. 27(b) which 
show the spectra for several time segments all 350 microsecs long but starting at different times, 
do not exhibit this effect, possibly because there are two pulses within the selected time segments, 
as shown by the magnitude spectra given in Fig. 28 for the same segments. Above 70 microsecs, 
the low frequencies dominate the spectrum, but both high and low components are present. 

By contrast, the phase spectra for a fixed time segment from a few different sensors (at dif- 
ferent distances from the impact) are shown in Fig. 29 (a) and (b) for Tests 04/06-3 and 04/14-3, 
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0 0.2 0.4 0.6 0.8 1 


AlfCULU fKIQUINCT (BADS/MICROSEC) 

(b) Main Pulse (150 ■ 250 microsac) 

Figure 25. Comparison of Magnitude Spectra 
for First and Main Pulses 
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(a) First Puis* (50 • 150 microsac) 



(b) Main Pulss (150 • 250 microsec) 


Figure 26. 


Comparison of Phase Spectra for First and Main Pulses 
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ANGULAR FREQUENCY (RAD8/MUSEC) 

(a) Phase calculated by Principal Value* of Arc Tangent 



ANGULAR FREQUENCY (RADS/MUSEC) 

(b) Continous Phase Calculated by Eliminating Jumps 
Note: Frequency scales are different in (a) and (b) 


Figure 27. Illustration of Continous Phase Calculation for 

Several Time Segments (Test 03/06, Sensor No. 4) 
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0 0.4 0.8 1.2 

ANGULAR PREQUERCT (RADS/MUSEC) 


Figure 28. Magnitude Spectra of Several Time Segments (Test 03AJ6, Sensor No. 4) 
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respectively. The slopes of these spectra are related to the distance of the sensors from impact, as 
shown in Fig 30. The frequency derivative of phase is determined by the group velocity, since 

d0/dco = (dk/dco) x - to = (1/Cg) x - to 

From the figure the data fit moderately well to a line drawn with a group velocity of 5.2 
mm/microsec, which is the velocity of the first wave. Since this is essentially non-dispersive, the 
group and phase velocities are the same. 

The magnitude spectra for two different impactors, a steel BB (Test 04/06-3) and a nylon 
sphere (Test 04/14-1) are shown in Fig. 31. Sensors no. 1 in (a) and no. 2 in (b) are at about the 
same distance of 540 and 570 mm. These two spectra show that the nylon produces higher fre- 
quencies than the steel. This is expected on the basis of elastic impact theory, since the nylon 
sphere would bounce off early, after only one wave reflection, being of lower accoustic impedance 
then the aluminum, whereas the steel BB would remain in contact for a long time. The other 
spectra in the two tests somewhat confirm this. 

5.4 Feature Analysis 

Another step in analysing the data was the use of Feature/Cluster Analysis for pattern recog- 
nition. The intent was to evaluate the potential use of the method for discriminating types of im- 
pact, and possibly for quantitative measures such as impact speed. A software package for ex- 
tracting features from waveform data, the 4060 System (reference 8) developed by General Re- 
search Corp. for their ALN 4060 Flaw Discriminator, was used to analyze some of the wave- 
forms. 

The software produces numerical values of 89 quantities such as peak amplitudes, total en- 
ergies, times between peaks, spectral maxima, and statistics of peaks. In evaluating its useful- 
ness for this application, the definition of the features extracted from an input signal have been re- 
viewed, and collected into groups considered appropriate, as given in Table 5. 

Some features, such as times of certain events, should be dependent on sensor position and 
some should not. Others should depend on impactor characteristics including speed or on plate 
response such as perforation. The features for sensors at three distances in Test 2/10 were re- 
viewed briefly for correlation to sensor position, and 6 have been identified as plotted in Figure 
32. These are described as: 

No. 37 (ALN FI) - 25-90% rise time No. 64 (ALN F34) - Power in first 1/8 of power spectrum 
No. 54 (ALN F13) - Phase coherence No. 60 (ALN F42) - Power in first 1/4 of power spectrum 
No. 55 (ALN FI 4) - Phase coherence 
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0 0.2 0.* 0.0 0 °- 2 °- 8 
ANGULAR FREQUENCT (RADS/MUSBC) ANGULAR FREQUENCT (RAD3/MUSEC) 


(a) Test of 04/06-3 (With Steel Impactor) <b) Twt of 04/1 4-3 (With Nylon Impactor) 

Figure 29. Phase Spectra for Several Sensors at Different Distances in Two Tests 



Distance, mm 


Figure 30. Slopes of Phase Curves for Two 
Tests as a Function of Distance 
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(a) Spectrum for Steal Impactor (Test of 04/06-3) 



(b) Spectrum for Nylon Impactor (Test of 04/14-3) 


Note: 

Spectra shifted vertically 
to separate curves 


Figure 31. Comparison of Spectra for Two Different 
Impactors— Steel and Nylon 
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Table 5. Definition of ALN Features and their Values for 2 Sensors in ALN Features for 04/06-3 No. 2 and 3 


Name 

ALN No 

Ranges 

Inputs 

Test No. 
4063/2 

Sensor No 
4063/3 

1 Pk amp 

57 

0-2048 

9507 

6168 

2 Pk amp in W3 

78 

(0-200) 

9507 

4426 

3 Pk amp in W4 

79 

(507-2048 

973.6 

895.5 

4 Mean Amp 

59 

0-2048 

1383 

994.9 

5 “ “ in W1 

68 

0-210 

319.1 

84.89 

6 “ “ “ W2 

72 

190-2048 

1500 

1087 

7 “ “ “ W4 

82 

(507-2048 

6364 

2672 

8 Std Devn 

60 


1075 

789.5 

9 “ “ in W1 

69 


386.2 

11.28 

10 “ “ “ W2 

73 


1060 

770.6 

11 “ “ “ W4 

83 


1330 

BIsEiKFICflS 

12 Ratio of 2 pks 

61 


0.6693 

0.717 

13 “ ““ “ in W1 

71 


0.1044 

0.5153 

14 “ “ “ “ in W2 

75 


0.599 

0.4178 

15 “ ““ “ in W3 

81 


4 

61 

16 “ ““ “ in W4 

85 


84 

29 

17 No pks above thld 

67 

30 

121 

125 

18 “ “ “ frac 1 

65 

10% 

100 

103 

19 “ “ “ frac 1 in W1 

70 

it 

1 

0 

20 “ “ “ frac 1 in W2 

74 

it 

100 

103 

21 “ “ “ frac 1 in W3 

80 

it 

1845 

720.7 

22 “ “ “ frac 1 in W4 

84 

it 

819.9 

581.6 

23 “ “ “ frac 2 

66 

35% 

14 

22 

24 A1/(A1+A2) in entire wfm 

62 


0.599 

0.4178 

25 Awl/(Awl+aw2) 

77 


0.1921 

0.0173 

26 Aw3/(Aw3+aw4) 

89 


0.599 

0.624 

27 Env time integ 

26 


1920000 

1440000 

28 Avg bef pk * avg aft pk 

32 


0.338 

0.433 

29 ditto shifted 

33 


12 

16 

30 Time at main pk 

58 


235 

1495 

31 Time betw 2 large pks 

63 


1176 

1175 

32 Time betw pk in W1 & W2 

76 


29 

1447 

33 Time betw pk in W1 & W3 

86 


0 

1175 

34 Time betw pk in W1 & W4 

87 


1176 

462 

35 Time betw pk in W3 & W4 

88 


1176 

1637 

36 tl/(tl+t2) 

64 


0.1428 

0.8237 

37 Rise time 

1 


2 

7 

38 Fall time 

4 


13 

23 

39 Pulse width @ .25 of pk 

3 


18 

34 

40 “ “ “ .6 of pk 

2 


6 

14 

41 Centroid of entire wfm 

27 


1041 

1115 

42 2nd mom (Std Dev) 

28 


587.6 

540.2 

43 3rd mom (Skewness) 

29 


0.0318 

-0.0933 

44 4th mom (Kurtosis) 

30 


1.742 

2.06 

45 Effec pk curv 

31 


0.354 

0.451 
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Table 5. Definition of ALN Features and their Values for 2 Sensors in ALN Features for 04/06-3 No. 2and 3 
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Sensor position, ins 

Figure 32. Correlating Features of Waveforms 
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6.0 CONCLUSIONS 

The potential of a space debris measuring system, based on collecting panel response data 
from particle impacts, has been evaluated from two major points of view: the statistics of the space 
debris environment were analyzed to determine the probability of success of the system, and ex- 
perimental data from impact tests were used to determine the usefulness of the measurements of 
the system. 

Particle flux data for space debris and data on its velocity distribution were combined to pro- 
duce probability data for momentum of impacts. These were then applied to the panels of the 
Space Station Logistics Module showing that for a typical impact detection sensitivity of 0.1 
g.m/sec, the probability of an event occuring in the six months that the module is in orbit is about 
0.3 for each panel of about 1.4 sq.m. For all the panels on the entire module the measurement 
system will then produce a marginal amount of data. Increased sensitivity of the detector system is 
certainly possible, and no attempt has been made to improve it. New sensor susytems, such as 
embedded acoustic waveguides currently being developed offer such improvement. 

The waveforms from impact tests have been analyzed to evaluate two factors: the potential 
accuracy of locating the impact point and the feasibility of extracting impact characteristics from the 
waveforms. 

The accuracy of measuring arrival times is within about 3 microsecs., which represents 15 
mm of propagation at the wavespeed (about 5mm/microsec) of the first in-plane wave motions. 
This leads to an accuracy of location of about 75 mm through the use of the location algorithm 
using 4 waveforms. Repeatability from test to test and between tests is within the same tolerance. 

The response data is shown, by analysis of the spectrum of the response, to be separable 
into the two modes of plate response: in-plane and transverse (bending) motions. The spectral 
magnitude is concentrated at high frequency for the first waves and at low frequency for the bend- 
ing. 

The in-plane response is the faster and arrives first, and at almost constant wavespeed. It is 
therefore well suited for use in the location algorithm. Methods for automatically identifying the 
arrival of this motion were developed. 

The bending motion is slower, and highly dispersive: the wavespeed increases with frequen- 
cy so that the high frequency components arrive first followed by low frequency motions later, 
stretching the pulse in time. The spectral behavior of the phase of these motions was shown to be 
predicted by plate bending theory and to be related to the propagation distance. This data could be 
used in the location algorithm or in a procedure for determining the propagation wavespeed. 
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Analysis of two waveforms taken at approximately the same distance from the impact in two 
tests, one with a steel BB and one with a nylon sphere, showed that the latter had a higher fre- 
quency spectrum as expected from impact rebound theory. 

Preliminary use of feature analysis showed that waveforms could be distinguished on the ba- 
sis of statistical analysis when distinct characteristics existed, such as the distance of the sensors. 
The data set, developed initially for evaluating location procedures, was inadequately varied with 
regard to impact parameters such as mass, velocity, or momentum to evaluate the ability to deter- 
mine these quantities. 

The conclusions reached are that the use of sensors to measure the motions of a panel after 
impact can produce data on the location and characteristics of the impact, and that the Logistics 
Module panels provide a possible application for detecting the low probability events of concern. 
The Method requires further development of data processing, and, to extend the range of detection 
requires the use of more sensitive sensors. 
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Appendix A 
Waveforms for All Tests 
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Figure A-2 . Tests of 02/10/87: Two Repeats at 5 in, 7 in, and 10 in 
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Second Test 

Figure AS Tests of 03/20/87: Two Repeats With Two Sensors at 10 in 
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Four Transducer) at 12 in 


Figure A-7. Second Test of 03/25/87: Four Transducers at 45° Separations on 12 in 
Circle and Trigger Transducer at 0 in 
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Figure A- 10. First Test of 04/06/87: 
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, Figure A - 14. First Two Test of 04/14/87 
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Figure A-24. Third Test of 04/23/87 
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! Figure A-28. First Test of 05/04/87 
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Figure A-30. First Test of 05/05/87 
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APPENDIX B 

Space Debris Momentum Distribution 
Fritz Scholz 

Boeing Computer Services 

Introduction: In order to obtain more accurate information about the 
space debris that forms a hazard for the space station it is proposed to 
use the cylindrical logistics module of the space station to gather space 
debris impact data. By equipping various outside panels of the cylinder 
with sensors it is possible to measure time, location and the surface normal 
momentum component M of each impact. In order to calibrate the impact 
sensors appropriately it is useful to have some idea of the distribution of 
the momentum M, i.e. how likely is it to get a momentum outside a certain 
range. Such an assessment is possible only if certain probabilistic modeling 
assumptions are made. These assumptions are not unreasonable but repre- 
sent some simplification of reality. Further, once the model is assumed and 
the mathematical form of the momentum distribution function is derived, 
its evaluation still depends on the angle and diameter distributions of the 
debris particles. These distributions are known only imperfectly and it is 
the objective of this data gathering experiment to improve on that knowl- 
edge or confirm it. The next section discusses the probability modeling 
assumptions that underlie the derivation for the momentum distribution 
function. The formula for the distribution function is stated and a proof 
is given. The final section gives a formula for the expected number of hits 
of diameter d or greater on a given panel and for the full cylinder mantle 
some numerical values are given on a yearly basis for several values of d. 

The Probability Model: Each debris particle can be characterized by its 
diameter D and the angle 0 that its relative velocity vector forms with the 
space station flight direction. This relative velocity is then v(0) = 2v 0 cos(0) 
where i>o is the common orbital velocity of space station and debris particle. 

It is assumed that D and 0 are independent random variables with re- 
spective densities h(d) and t/;(0). The distribution function of D is then 
H(d) = P(D < d) = / 0 d h(x) dx and 1 — H(d) = P{D > d) is denoted by 
H{d). It is assumed that the density tv(0) of 0 is explicitly given. H(d) is 
given only indirectly through the debris flux measure F Deh [d ) which repre- 
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sents the average number of debris particles of diameter > d that will pass 
through a cylinder with cross section 1 in unit time (the cylinder standing 
perpendicular to the orbital plane). This flux definition is bigger by a factor 
of 4 than the tumbling plate flux, see [l]. On page 13 of [l] it was shown 
that 

XWMO) = (°) |cos(0 - «] + 

which implies that H(d) = Fp^d)/ F^^iO). There is a problem in that 
F Deb (d) is not given for d < 10“ 2 cm and extrapolation suggests that 
Fotb (0) = oo. A way out of this dilemma is to define a threshold e > 0 
such that diameters < € are no longer considered to be space debris. Then 
define H(d) = 1 — Fr)tb{d) / Fd^c) for d > e and H(d) =0 for d < e. What 
constitutes a reasonable threshold e is left open at this point. 

Since the sensor will measure the surface normal momentum component 
of the debris particle it becomes necessary to introduce a surface area ele- 
ment A of the cylinder. It is most convenient to consider surface elements 
that represent narrow vertical rectangular strips on the outside cylinder 
wall. Such a segment is characterized by the angle <f> that the ray from the 
cylinder center to strip centerline forms with the space station flight direc- 
tion. <f> may vary from 0 to 2tt. If r represents the cylinder base radius, h 
the height of the rectangular strip and rd<j> its base, then the area |A] of A 
is |A| « hrd<f>. To emphasize the dependence of A on <f> the notation A$ is 
used. The normal momentum component of a debris particle of diameter D 
that impacts A $ from direction 0 is then M = cD z v(9) [cos(0 — <f>)\ + where c 
is the proportionality factor that relates volume to mass and [cos(0 — <^)] + 
accounts for the surface normal part of the momentum which becomes zero 
whenever cos (0 — <t>) < 0. 

The random process of arrivals of debris particles with their respective 
random angles and diameters D is modeled by a marked Poisson process 
Nj, with intensity 

\f(t,0,d) = Jz? e i(e)|A*|[cos(0 - ^)] + ttf(0)h(d)J[o iO o)(t) 

which is defined on S = [0, oo) x [ — zr/2, tt/2] x [0, oo). For any subset C C S 
N${C) represents the random number of debris particles which hit A^ and 
which have characteristics ( t,0,d ) 6 C. Here t refers to the arrival t im e 
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of the particles. Further, for any C C S the random number i\^(C) has a 
Poisson distribution with mean 

H*{C)=[ A {t,0,d) d(t,$,d) 

jc 

and for any k = 1,2,3,... and disjoint Ci, . . . , C* C S the random vari- 
ables JV,(Cx), ..., N$(C)c) are independent. For more details on marked 
Poisson processes see [2], [3]. 

A given vertical rectangular panel B on the outside of the cylinder can 
be thought of as approximated by many. narrow strips A*, <j>i < 4> < fa, 
where <f > x and <f> 2 delimit the base arc of the panel B. For each narrow strip 
a marked Poisson process N, $ is modeled and it is further assumed that 
these Poisson processes for different strips are independent. This construc- 
tion seems reasonable enough since those particles that hit one strip on 
the outside are separate from the particles that hit any other strip on the 
outside. 

The Momentum Distribution Function: Having made these modeling 
assumptions it is possible to give the following formula for the distribution 
function of M, namely 

G(m) = P(M < m) = V/R , 

where 

y = C G |cos(s ■ '* )1+u,(i ’ ) * ) dei * 

and 

R = [ [ [cos(0 — <j>)] + w(0)d9d(f > . 

i J-r/2 

The proof for this will be given only for the momentum Mi of the first 
particle that hits panel B. It is the same for any particle but the derivation 
is somewhat messier in notation. 

Proof: If Mi is the surface normal momentum component of the first 
debris particle to hit panel B and Ti its time of arrival then 

G(m) = P(M\ < m) = lim P(Mi < m, rj < T) 

T — ►OO 
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and 

P{M X < m, Tl < T) = f; P{M X <m,T 1 € /,) 

»=1 

where is the interval [T(t — l)/n,Tt'/n) for t = 1, 2, n. Let 

° im = : 4 G Ti ' d ~ (ct;(0)[cos(0 -£)]+) ' 6 e J *} 

where = {6 € f— 7r/2,jr/2] : cos (6 — <j>) > 0}. C{ m delineates the 
characteristics of any particle arriving on A$ during / ,• with surface normal 
momentum component < m. Similarly, let 


C,-i = {(f,M) : t e [0,T(t - l)/n), d > 0, 9 e J+} 

Then P[M X < m, t x £ /,) can also be written as 

P{N 4> {C i - 1 ) = 0 V<£ e [^ 1 ,^ 2 ] > 1 for some ip € [fa, fa]) 

= P(i\fy(C,-_ 1 ) = 0 V<£ € [^i,fo])P(JV^(C {m ) > 1 for some ip e [fa, fa]) 
= II P{N*(Ci- 1 ) = 0 )P(N i> {C im ) > 1 for some rjj € [fa, fa]) 


= e xp(-J3 ^{ c i-i)) p i N A C im) > 1 for some ip e [fa, fa]) 

4 


where 


t — 1 rr\2 

Vt(Ci-i) = ~T / P, ei (e)K|[cos(0 - *)]+«,(*) d6 

Tl J-r/2 

so that (j,$(Ci-i) can be approximated by 

1 —^TF deb (£)rh [*' r /2 [cos(0 - <f>)] + w{6) d0d<f> = ^ TF deb (e)rhR . 

n J-t /2 n 

Since 


P(^(a-m) > 1 ) = P{Ni,{C im ) = 1) -f P(iV4C im ) > 1) 
= MCim) exp(-^(C im )) + P(N^(C im ) > 1) 
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and 


P(N 4 C im ) > l) — 1 — exp( fii/)[Cir, t)) < n^(Ci m )/2 

and since 


= - r f Drt ( e )|X*|[cos(«-^)] + u.( 9 )F 
M J-ic/2 


m 


l/3> 


v ct>(0)[cos(0 — V’)] + y 




which is of order l/n, we have 

P(JV*(C fm ) > 1 ) = P(JV*(C im ) = 1 ) + 0 ( 1 /™*) 

and 

P(W,(C im ) = 1 ) = /x,(C fm ) exp (-^(C im )) = 0 (l/n) . 

Thus (using Bonferroni’s inequality) we have 

P(N+(C im ) > 1 for some ijt € [<f> i,&]) = 2 Z P(^(C im ) = 1 ) + 0 (l/n 2 ) 


and 


J3P(JV*(C im ) = x ) = H^(C'«m)exp(-^(C, rn )) 

V> ^ 


—rhF Dtb (e) [ [ / [cos( 0 -^)] + u>( 0 ).ff 

n J-*i 2 


m 


l/3\ 


cv(0)[cos(0 — 


dddip 


= ~rhF Dei {e)V . 

n 


Hence 


P(M i <m, Ti < T) = exp 

i=l 


* 1 -rJ , j,.»(e)rfcB| -rhF I)e6 (€)7+nO(l/n 2 ) 

/ n 


n 


which, as a Riemann sum, approaches 

/ exp(-xTF Dei (e)rhR)dxTrhF D ' b (e)V = j^(l-ex-p(-TF Deb (e)rhR)) 
as T — ► oo. Hence P(M X < m) = V/R as claimed. 
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Expected Number of Impacts: In order to gauge the feasibility of the 
experiment it is necessary to find out how many impacts of space debris 
particles of size d or greater can be expected for a given panel B. To 
obtain this expected value one proceeds as before by viewing the panel as 
approximated by many narrow panel strips and by adding the expected 
values for each of those. To get the expected value for a panel strip let 


C T = {{t,d,d) : d>d 0 ,t<T,eeJt}, 

then 

hAC t ) = TH(d 0 )rhd<f>F Deb (e) f [cos(0 - 4>)] + w(9)d9 

= Trhd<f>F Dei (do ) / [cos(0 - <f>)] + w{9)d9 
J-x/2 

represents the expected number of debris impacts of diameter > d 0 during 
a time period of length T on the strip corresponding 4>. By summing this 
over <j> € [<^ 1 ,^ 2 ] one obtains the corresponding expected number for the 
full panel B as 

TrhF Dei (do) f I*' 2 [cos(0 - <j>)} + w(9)d9d<f> , 

H x J~ir/2 

which for <f>i = 0 and <f >2 = 27 r simplifies after interchange of integration 
to 2Tr/iFu e j(<io), which fits in with the used definition of Fu t j since 2rh 
represents the cross section of the cylinder. Note that the expected number 
of impacts is independent of the threshold e provided do > e. 

Taking the flux data for 400 km altitude as given in [4] for a tumbling 
plate and noting the factor 4 increase of Fd^ over this tumbling plate flux 
one can read off the following formula for Fceb- 

F D eb{d) = 6.34 • 10 -6 • d~ 2A for d < 1. 

Since the cross section area of the logistics module is about 23.7 m 2 one 
can expect about 9.5, .038, and .00015 yearly impacts of diameter > .01, 
> .1 and > 1 respectively, on the outside mantle of the logistics module. 

Also of interest is the expected number of impacts with surface normal 
momentum component > m. This expected value is obtained as before by 


c -3L 
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integrating the expected values over narrow strips. The expected value for 
a narrow strip is obtained by letting 

c T , m = ■ t<T,d> ( cu(0 ^ cos ^_^| + ) ’ J *j 

and by computing 

H[CT„)=TF D 4')J^ lt (( CT( , )[co ^_fl |+ 

Summing this over <f> € [$ 1 ,^ 2 ] one obtains for the expected number of im- 
pacts of momentum components > m on the panel B the following formula: 

TrkKM £ 7 * [«.(« - «♦-«>* ded * 

which reduces to 

w: c, ; ) iH * ■ 

Note that the dependence on e has disappeared as long as e is chosen smaller 
than any possible argument in the F^eb function. This depends mainly on 
the choice of m. Since Fr>cb{d) is not given for d < 10 _2 cm one still has 
to make the decision on how to extrapolate Fz> e t,(d) below 10 -2 cm. One 
approach is to take the linear extrapolation on the log-log scale which leads 
to 

F Deh {d) = 6.34 • 10 -6 • d~ 2A for d < 1. 

The other approach is to truncate Fz> t b{d ) for d < 10 -2 cm by setting 
F Dtb {d) = F£> < j(10 -2 cm) for d < 10 ~ 2 cm. The rationale behind that is 
to treat all particles with diameter < 10 -2 cm no longer as debris. 
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